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THE VIEWING ANGLES OF BROAD ABSORPTION LINE VERSUS UNABSORBED QUASARS 
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ABSTRACT 

It was recently shown that there is a significant difference in the radio spectral index distributions 
of broad absorption line (BAL) quasars and unabsorbed quasars, with an overabundance of BAL 
quasars with steeper radio spectra. This result suggests that source orientation does play into the 
presence or absence of BAL features. In this paper we provide more quantitative analysis of this 
result based on Monte-Carlo simulations. While the relationship between viewing angle and spectral 
index does indeed contain a lot of scatter, the spectral index distributions are different enough to 
overcome that intrinsic variation. Utilizing two different models of the relationship between spectral 
index and viewing angle, the simulations indicate that the difference in spectral index distributions 
can be explained by allowing BAL quasar viewing angles to extend about 10° farther from the radio 
jet axis than non-BAL sources, though both can be seen at small angles. These results show that 
orientation cannot be the only factor determining whether BAL features are present, but it does play 
a role. 

Subject headings: quasars: general, quasars: absorption lines 



1. INTRODUCTION 

A long-time popular explanation for the presence of 
broad absorption lines (BALs) seen in approximately 
20% of quasar spectra (Knigge et al. 2008) has been 
a simple orientation model, in which BAL quasars are 
seen only from a more "edge-on" perspective, or at larger 
viewing angles (Elvis 2000). Understanding the geome- 
try of the outflows producing these lines is an important 
part of modeling the role they play in the evolution of 
the quasar itself, as well as the effects they may have on 
the surrounding environment and host galaxy via feed- 
back effects. For example, it has been shown that it is 
possible for AGN feedback to affect star formation rates 
in the host galaxy, both theoretically (Hopkins & Elvis 

2010) and now it seems observationally (Cano-Diaz et al. 
2012). 

The similarity of the emission lines in BAL and non- 
BAL quasars (Weymann et al. 1991), as well as their 
optical polarization properties (Ogle et al. 1999) have 
been used to support simple orientation models. How- 
ever, this scheme fails to explain various other observa- 
tions, particularly at radio frequencies. Short timescale 
radio variability has been identified in around 20 BAL 
quasars, which is argued to indicate a viewing angle near 
the radio jet axis to explain the derived brightness tem- 
peratures (Ghosh & Punsly 2007, Zhou et al. 2006). In 
general, BAL quasars are more compact than non-BALs 
in radio maps (Becker et al. 2000), at least at low to 
intermediate resolution, and in small samples they do 
not show a significant difference in radio spectral index 
distribution compared to non-BAL sources (Becker et 
al. 2000, Montenegro-Montes et al. 2008, Fine et al. 

2011) . Radio spectral index (a; S v oc z/ a , where S u is 
the radio flux and v is the frequency) is generally con- 
sidered an orientation indicator, with steeper spectrum 
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(a < —0.5) sources seen more edge-on because they are 
dominated by optically thin lobe emission, which has a 
steep spectrum due to little synchrotron self-absorption. 
More pole-on sources are dominated by core emission 
because of relativistic beaming effects, and have a flatter 
spectrum because they are optically thick and thus sig- 
nificantly self-absorbed. Due in large part to these obser- 
vations, other explanations based on pure evolution (e.g. 
Gregg et al. 2006), and not orientation, have also come 
into favor. 

DiPompeo et al. (2011) expanded greatly the num- 
ber of BAL quasars with multi-frequency radio data and 
presented a sample of 74 BAL quasars, along with a sam- 
ple of 74 individually matched unabsorbed quasars, with 
flux measurements at 4.9 and 8.4 GHz (observed frame) 
from the Very Large Array (VLA) /Expanded Very Large 
Array (EVLA). These data provided quasi-simultaneous 
flux measurements to remove the effects of radio vari- 
ability in the spectral index measurements. The distri- 
butions of | were significantly different, with BAL 
quasars showing an overabundance of steep spectrum 
sources, but both samples show a wide range of spectral 
indices. Analysis of other measurements of a (includ- 
ing a fa-, a simple linear fit to available literature fluxes 
at various frequencies, in addition to the new measure- 
ments), a variety of statistical tests, and a restriction to 
only unresolved sources in both samples all also show 
that the difference is present, although the significance 
does vary. Some of this variation could be due to the fact 
that the other measures of a included non-simultaneous 
flux measurements, or due to variation in the number of 
sources included in the tests. Because the two samples 
are one-to-one matched in redshift, among other proper- 
ties, use of rest-frame spectral indices will not effect the 
results. 

These results indicate that while BAL quasars likely 
span a range of orientations, viewing angle does plays a 
role in their presence. The next step, presented here, is to 
quantify this difference and provide the most likely view- 
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ing angles to these sources, at least in a general sense. 
Our aim is to test if the difference seen in spectral index 
distributions can be explained by differences in viewing 
angle, and we recognize that it may be possible to de- 
velop more sophisticated models and simulations in the 
future. 

2. THE a-9 RELATIONSHIP 

We have based our modeling off of two relationships 
between a and viewing angle (9, defined as 0° along 
the radio jet axis); one purely observational, and one 
from semi-empirical simulations. We began developing 
the observational relationship with the sample of Wills 
& Brotherton (1995), which included the 29 quasars of 
Ghisellini et al. (1993) with viewing angles calculated 
from superluminal motion seen in VLBI maps, as well 
as 4 additional sources from Vermeulen & Cohen (1994) 
with superluminal motion measurements and X-ray data 
at 1 keV available. We then searched NEE0 for multi- 
frequency radio fluxes in order to measure a for these 
sources. We initially attempted to use fluxes at 4.85 and 
8.4 GHz to better match the observations in DiPompeo 
et al. (2011), as well as a linear fit to all available ra- 
dio fluxes, but given the small number of sources and 
the scatter in the a-9 relationship we were unable to get 
a reasonable fit using these data. One complication is 
that in general the radio flux measurements are not si- 
multaneous, and variability almost certainly exaggerates 
the scatter since the majority of these sources are seen 
from small viewing angles. In the end we settled on using 
the spectral index between 15 and 8.4 GHz, which was 
available for 27 of the 33 sources (26 from Ghisellini et al. 
(1993) and 1 from Vermeulen & Cohen 1994), in order to 
build a useable model. We note that a necessary assump- 
tion in this analysis is that on average the radio spectra 
are reasonably approximated by a power law over a large 
frequency range, which may be an oversimplification in 
some cases. However, if you inspect the spectra presented 
in DiPompeo et al. (2011) for both BAL and non-BAL 
samples, when data is available at more frequencies the 
spectra are often well behaved across a large frequency 
range. 

We adopt the values of the viewing angles from Ghis- 
ellini et al. (1993) and follow their method to calculate 
the viewing angle for the additional Vermeulen & Co- 
hen (1994) source. Two sources that were large outliers 
(1830+285 and 1845+797, from Ghisellini et al. (1993), 
which both have highly inverted radio spectra) were ex- 
cluded in order to allow us to make a fit to the data that 
could possibly reproduce our observed range of spectral 
indices. The final list of sources, their viewing angles and 
their spectral indices are given in Table [TJ 

While the relationship is likely more complex in reality, 
we made a simple linear fit to the data. Above a view- 
ing angle of 30° we assume a flat relationship where all 
the scatter is due to intrinsic variation in radio sources 
for several reasons. First, there is no observational data 
above 30° to constrain the shape of the relationship. Sec- 
ond, if we simply extend the linear relationship, we would 
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see large numbers of sources with extremely steep spectra 
(a extending to —3 or —4), which is simply not observed 
in most radio sources. Finally, although it is dependent 
on the Lorentz factor of the the emitting material, the 
increase in observed flux due to Doppler boosting is ex- 
pected to be small at larger viewing angles. Because vari- 
ability almost certainly exaggerates the scatter in the re- 
lationship, we use the standard deviation of the spectral 
index distribution of the quasars in the 3CR catalog (as 
presented in Smith & Spinrad 1980) to model the scat- 
ter. Because this catalog was built using low frequency 
data, it is dominated by lobe emission and therefore con- 
sists of mostly steep spectrum, edge-on sources. Shown 
in the top of Figure Q] is the data used from Ghisellini 
et al (1993) and Vermeulen & Cohen (1994) and the fit 
used to model the a-9 relationship; the bottom of Fig- 
ure Q] shows the result of our Monte-Carlo simulation of 
10,000 randomly distributed jet viewing angles and their 
spectral indices based on the above model. 

The semi-empirical model used is from the simulations 
of Wilman et. al (2008), which is the extragalactic por- 
tion of the SKA Simulated Skies (S 3 ) project. We pulled 
from their simulations the fluxes at all available frequen- 
cies (151 MHz, 610 MHz, 1.4 GHz, 4.86 GHz, and 18 
GHz, all observed frame) for all FR I and FR II type 
sources between redshifts of 1.5 and 3.5 and with values 
of Si. 4 > 10 mjy Ideally we would use the spectral index 
between 4.86 and 18 GHz from these simulations; how- 
ever, there are issues with their high frequency results. 
The source populations are drawn from a 151 MHz lu- 
minosity function and extrapolated to high frequencies, 
which almost certainly does not accurately predict the 
high frequency source population (see for example, Ma- 
hony et al. 2011). Additionally, there may be problems 
with the functional form assumed for the SEDs of radio 
cores, causing them to fall off too steeply at high fre- 
quency (Wilman, private communication). This causes 
the maximum spectral index found to drop well below 
as you go to higher frequencies, which has no bearing 
on physical reality. The simulations also apply a lower 
limit on the lobe spectral indices at a = —0.75, which is 
problematic because we clearly see a high fraction of our 
sources with steeper spectra than this (the steep spec- 
trum overabundance in BAL sources is most prominent 
between —2 < a < —1). However, if we again sim- 
ply assume that in general radio spectra of these objects 
obey a simple power law, we can use the lower frequency 
spectral index distribution from these simulations and as- 
sume that a similar relationship holds between 9 and a 
regardless of which part of the radio spectrum is consid- 
ered. Given the scatter in the relationship it is unlikely 
that this assumption will significantly effect our results. 

The relationship between 9 and ctgio from these sim- 
ulations is shown in the upper panel of Figure [2j In 
order to use this in our simulations, we made two fits 
to different ranges of 9. First, we did a simple linear fit 
from 10 < 9 < 20, where the data can reasonably be 
approximated as linear. This allows us to have sources 
with steeper spectra than a — —0.75, which is needed to 
reproduce what is actually observed despite the limit ap- 
plied in the simulations. The distribution in a is then cal- 
culated by normalizing the data by the fit so all the vari- 
ation is in a, assuming the distribution is Gaussian and 
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calculating the standard deviation. The data is clearly 
non-linear at small 8, and so between 5 < 8 < 20 we also 
make a separate linear fit in log(a + 0.75)-$ space. Our 
final model is then based on a combination of these fits; 
below viewing angles of 10° the fit in logarithmic space is 
used, and above 10° the linear fit is used. Both lines are 
shown in the upper panel of Figure The lower panel 
of Figure [2] shows the result of our Monte-Carlo simu- 
lation of 10,000 randomly distributed jet viewing angles 
and their spectral indices based on this semi-empirical 
model. 

The overall general shape of the two models used is 
similar; however, the amount of scatter in the observa- 
tional model is significantly more than that in the empir- 
ical model. We also do not consider the possible effects 
of redshift in either model, though it is possible that 
the wide range of redshifts in the observed sample could 
have an effect on the spectral index distributions. Again 
however, since the BAL and non-BAL samples are well 
matched in redshift, we feel that any effects of redshift 
will not change the general trends found. 

3. SIMULATIONS 

The simulations were performed using IDL. We first 
generate a random vector in 3-D space, utilizing IDL's 
uniformly distributed random number generator to cre- 
ate x, y, and z coordinates between -1 and 1, which are 
combined to create vectors within a "unit cube" . This 
will produce a distribution of vector directions biased to- 
ward the corners of the cube, and so we then reject any 
vector that does not fall within the unit sphere (has a 
magnitude greater than 1), resulting in vectors uniformly 
distributed in random directions in 3-D space. This vec- 
tor is taken to represent one side of a bipolar radio jet. 

We next choose a viewing direction along the z-axis, 
and compute the viewing angle to the randomly gener- 
ated vector. If the vector is pointing in the negative z 
direction, we reflect it about the origin by multiplying all 
of its components by —1, to represent the opposite side 
of the bipolar jet, so that viewing angles range from 
to 90°. Once the viewing angle is known, we can assign 
a value of a based on one of the models in the previous 
section. This is done utilizing IDL's normally distributed 
random number generator, so that the assigned values of 
a are normally distributed about the model fit with the 
assumed standard deviation. The bottom panels of Fig- 
ures [D and El were generated by repeating this process 
10 4 times. 

In order to compare the Monte-Carlo simulations to 
the real data and quantify the viewing angles to BAL 
and normal quasars, we used the simulation to virtu- 
ally repeat our observational program while restricting 
what source viewing angles were allowed. We system- 
atically stepped through all possible ranges of allowed 
viewing angles between some 8 m in and 8 max , assigning 
values of a for each random viewing angle between the 
allowed 8 m i n and 8 max until there were 74 virtual "ob- 
servations" (to compare to the 74 objects really observed 
in each sample). We then compared the distribution of 
simulated spectral indices to the observed spectral index 
distributions for the BAL and non-BAL samples, using 
both K-S and R-S tests, to check whether the real data 
were well reproduced. Our criterion was a p- value of less 
than 0.05 to indicate that the distributions are from a 



different parent population, and are not well matched. 
We then repeat this experiment 10 5 times, in order to 
get a statistical sense of how often the real observations 
are well reproduced for a given range of viewing angles. 
Once this is complete, the allowed viewing angle range 
is changed and the process is repeated until all possible 
8 ranges are tested. We apply an upper limit to 8 max 
of 45° as above this value it is likely that most quasars 
are obscured by dust (Barthel 1989). As the results will 
show in the next section, this upper limit has no effect 
on our findings. 

4. RESULTS 

The probabilities of each set of simulations producing 
a match (based on the K-S tests) to the observed data 
as a function of 9 m i n and 8 max are shown in Figures IS} 
[6l The z-axis represents the percentage (out of the 10 5 
runs) that the resulting a distribution matched the ob- 
served distribution, the x-axis is 8 m i n and the y-axis is 
8 max- In each figure, the left panel shows the probabili- 
ties for the non-BAL sample, and the right panel shows 
the probabilities for the BAL sample. Figures showing 
the probability distributions using an R-S test for com- 
parison are not included to save space, since the results 
are not significantly different than using K-S tests. The 
results are all summarized in Tables [2] (K-S test com- 
parison) and [3] (R-S test comparison); Column (1) is 
the sample that the simulations were compared to, Col- 
umn (2) is the model used in the simulations ("Obs" is 
the model based on the sample in Wills & Brotherton 
(1995), "Emp" is the model based on the semi-empirical 
simulations of Wilman et al. (2008)), Column (3) is the 
maximum probability reached, Column (4) is the value 
of 8 m i n at P m ax, and Column (5) is the value of 8 max at 
Pmax- The final two columns, (6) and (7) to the right 
of the vertical line, indicate the values of 8 m in and 8 max 
with the largest separation that also have a probability 
of reproducing the observed results of greater than 90%. 

In general, the simulations are able to accurately re- 
produce our observations, with more than half of them 
having probability distributions peaking at or above 99%. 
The only exception is trying to reproduce the a/u dis- 
tribution for non-BALs using the semi-empirical model 
and comparing using a K-S test, where the probability 
peaks at only 71%. However, the shape of the probability 
distribution and the location of the peak is still basically 
the same as other simulations. The probability distribu- 
tions are also generally quite flat in the m in direction, 
except for the cases where the semi-empirical model and 
a K-S test are used; in those situations the distributions 
are well-peaked. This flatness is the reason we include 
the final two columns in Tables [2] and [3l as well as as the 
fact that the viewing angle ranges suggested by consid- 
ering only the peaks can often be unrealistically narrow. 
While there is some variation in the results depending 
on which model or statistical test is used, they are usu- 
ally quite consistent. In all cases, BAL quasars can have 
small viewing angles, but always extend farther from the 
radio jet axis when compared to non-BAL quasars. If 
we take an average of all simulations, we find an average 
viewing angle range of 0-22° for non-BALs and 1-32° for 
BAL quasars. 
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5. DISCUSSION 

With maybe the exception of the simulations using the 
observational model of the a-9 relationship and compar- 
ing to a fa with a K-S test, the general trend seen is that 
BAL quasars cover the same range of viewing angles as 
normal quasars, but extend about 10° farther from the 
jet axis. These results suggest that objects such as "po- 
lar" BALs (Ghosh & Punsly 2007, Zhou et al. 2006) are 
indeed real and that for small to intermediate viewing 
angles it is possible to observe a quasar either with or 
without BALs, but at the largest viewing angles one will 
only see BAL features. While constraining the a-9 rela- 
tionship can be difficult, it is interesting and encouraging 
that the results are quite similar using either of the two 
models or statistical comparisons here. We suggest that 
more effort be placed on constraining this relationship 
(both observationally and theoretically) in the future, as 
sample sizes grow to significant enough numbers to use 
radio spectral index as a statistical orientation indicator 
when no other method is available. 

One problem that is readily apparent in these results 
is that none of the samples extend to viewing angles as 
large as one might expect from, for example, the results 
of Barthel (1989). This could indicate a problem with the 
models used. However, the fact that the general result 
seems mostly independent of which model is used, we 
believe that while a change in the model may affect the 
absolute numbers it will not change the main conclusion. 
Also, as mentioned in DiPompeo et al. (2011) there may 
be a slight biasing of our original sample toward more 
face-on sources (due to the requirement that all sources 
have 1.4 GHz FIRST fluxes greater than 10 mJy), but 
this bias should effect both BAL and non-BAL quasars 
equally So while it is possible that 9 max may be higher 
for both samples, the end result should remain the same. 

Another consideration in this analysis is whether we 
can interpret the shape of the radio spectrum in a sim- 
ilar way in unresolved and extended radio sources. To 
build models of the a-9 relationship, we are required to 
utilize large scale, extended radio sources because these 
are the only sources in which 9 can be measured directly, 
at least in sufficient numbers to build a useable model. In 
contrast, the samples of DiPompeo et al. (2011) consist 
of high luminosity, generally compact radio sources (86% 
and 78% of BAL and non-BAL quasars, respectively, at 
5" resolutions), although there is generally not data at 
high enough resolution or enough frequencies to classify 
them as belonging to the special classes of compact steep- 
spectrum (CSS) or especially gigahertz peaked-spectrum 
(GPS) sources. It is possible that the shape of the spec- 
trum in the samples is affected by evolution, as is the 
case for CSS sources, and not just geometry. However, 
because CSS sources are selected to have steep spectra 
by definition, a comparison of spectral index distribu- 
tions to test this theory would not be useful. Instead, 
more radio data are needed, in particular at lower fre- 
quency, to look for turnovers in the radio spectra before 
comparisons with this class of object can be reasonably 
made. 

There are also indications that compact sources may 
have lower bulk velocities and thus lower Doppler fac- 
tors (e.g., Polatidis & Conway 2003). If this is the case 
in these quasars, it is possible then that the core compo- 



nent is not enhanced at the same level as other sources 
and therefore the lobe component remains dominant to 
smaller viewing angles. However, the fact that DiPom- 
peo et al. (2011) still see a significant difference in spec- 
tral index distributions even when restricting the BAL 
and non-BAL samples to only unresolved sources sug- 
gests that there are still significant orientation effects. 
Additionally, VLBI studies of BAL quasars (such as Doi 
et al. 2009, Kunert-Bajraszewska et al. 2010) do in fact 
show that compact BAL quasars have highly collimated 
radio jets, which means that viewing angle should effect 
the steepness of the radio spectrum in these sources in a 
similar way as more extended sources, as supported for 
example by the results of Jiang & Wang (2003). 

Regardless of the complications mentioned above, 
these simulations should lay to rest any claims that BAL 
quasars are only seen edge-on, and move the argument 
away from simple dichotomies toward a more complex ex- 
planation that includes both orientation and other fac- 
tors. Of course with the data used here, we can only 
make these claims for radio-loud BAL quasars, though 
the similarities in other geometrically dependent prop- 
erties between radio-loud and radio-quiet BAL quasars 
(for example, optical polarization; DiPompeo et al. 2010) 
suggest they may extend to radio-quiet BALs as well. It 
is clear that to explain the spectral index distributions 
seen, there needs to be a large overlap in viewing angles 
to BAL and non-BAL quasars. One simple explanation 
could be that BAL winds are launched at a variety of an- 
gles in different sources, or that they have a wide range 
of opening angles. It is also possible that BALs can be 
explained by a combination of orientation and evolution. 
The evolutionary schemes put forth for example by Gregg 
et al. (2002, 2006) suggest that BAL quasars begin en- 
shrouded in a cocoon of gas and dust, which is blown 
out as radio jets develop and luminosity increases. We 
may be able to draw on this description with these re- 
sults, by making the clearing out of the gas and dust 
orientation dependent. It is possible that the polar re- 
gions around the developing radio jets are cleared out 
first, which would explain the low numbers of seemingly 
"polar" BALs found as well as the spectral index distri- 
bution we observe. The most equatorial regions then are 
never completely cleared of absorbing material (or are 
replenished by a disk wind), causing only BAL quasars 
to be seen beyond a particular viewing angle. Whether 
a BAL is seen at small or intermediate angles is deter- 
mined by the evolutionary status of the quasar, but in 
the latest stages of the quasar lifetime the presence of 
BALs is only orientation dependent. We will attempt to 
develop this picture further in future papers, considering 
more factors than just the spectral index distributions 
modeled here. 

6. SUMMARY 

We have performed Monte-Carlo simulations of ran- 
domly oriented bi-polar radio jets and their correspond- 
ing radio spectral indices, in order to quantify the results 
of DiPompeo et al. (2011) that show a significant differ- 
ence in the spectral index distributions of BAL versus 
non-BAL quasars. The simulations utilize two different 
models of the relationship between spectral index and 
viewing angle; one based on observations and one based 
on semi-empirical simulations. By limiting the allowed 



source viewing angles and comparing the resulting spec- 
tral index distributions with our observations, we can 
constrain the viewing angles to BAL and unabsorbed 
sources. The results are mostly independent of which 
spectral index-viewing angle model is used, which spec- 
tral index measurement they are compared to, and which 
statistical test is used to compare the simulations to re- 
ality, at least in a general sense. We find that there is a 
large overlap between viewing angles to the two samples, 
with both probably extending all the way down to 0° or 
along the jet axis, which supports the claim that there 
are BAL quasars with polar outflows. However, viewing 
angles to BAL sources generally extend farther from the 
jet axis, with about a 10° span in which BALs will always 
be seen. 

It is clear that a simple orientation-only model for BAL 
quasars cannot explain all of their observed properties, 
though orientation does play a role. We need to move 
away from simple dichotomies in order to fully under- 
stand them, and it is likely that a combination of pre- 
vious models will prevail in explaining this important 
subclass of quasar. 

This research has made use of the NASA/IPAC Extra- 
galactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy, under contract with the National Aeronautics and 
Space Administration. M. DiPompeo would like to thank 
the European Southern Observatory for funding visits to 
collaborate with C. De Brcuck, as well as the Wyoming 
NASA Space Grant Consortium for providing funding 
to perform all of the radio observations. M. DiPompeo 
would also like to thank Dr. Ruben Gamboa of the Uni- 
versity of Wyoming Computer Science Department for 
assisting with the simulations, as well as Dr. Adam My- 
ers for providing computing facilities to run the simu- 
lations. Finally, we thank NRAO for their support in 
performing the radio observations that are the basis of 
these simulations; The National Radio Astronomy Ob- 
servatory is a facility of the National Science Founda- 
tion operated under cooperative agreement by Associ- 
ated Universities, Inc. 
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TABLE 1 

Sources Used in the Observational Model 



Source 


(°) 


„8.4 
"15 


Source 


(°) 


„8.4 
"15 


0016+731 


5.0 


0.91 


1040+123 


17.0 


-0.80 


0106+013 


3.0 


0.47 


1150+812 


13.2 


0.46 


0153+744 


26.1 


-1.53 


1156+295 


2.2 


1.12 


0212+735 


6.2 


0.28 


1226+023 


6.3 


-1.04 


0333+321 


3.0 


1.02 


1253-055 


3.2 


-0.06 


0430+052 


9.9 


0.37 


1641+399 


5.6 


0.56 


0552+398 


20.2 


-1.11 


1828+487 a 


6.8 


0.74 


0615+820 


21.7 


-0.56 


1928+738 


7.4 


0.11 


0836+710 


4.7 


0.21 


2134+004 


0.1 


-0.53 


0850+581 


12.5 


-0.93 


2223-052 


1.6 


-0.31 


0906+430 


0.9 


0.69 


2230+114 


4.0 


0.52 


0923+392 


4.4 


0.42 


2251+158 


5.8 


0.82 


1039+811 


24.6 


-0.32 









Note. — All sources and viewing angles arc taken from Ghiscllini et 
al. (1993) unless noted otherwise. The spectral index afg 4 is measured 
from data gathered via NED. 

a This source is from the sample of Vermeulen &; Cohen et al. (1994). 
The viewing angle was calculated using the method of Ghiscllini et al. 
(1993), using x-ray data from Wilkes et al. (1994) and additional radio 
data from Polatidis et al (1993). 



TABLE 2 

Simulation Results Using K-S Tests 



Sample 


Model 


a 


Pmax 


@min,P max t 


max, Pmax 


9min,P>».<3 


&max,P>0.9 


non-BAL 


Obs 


„4.9 
"8.4 


0.99 


12 


19 





24 


BAL 


Obs 


„4.9 
"8.4 


0.99 


13 


27 





34 


non-BAL 


Emp 


„4.9 
"8.4 


0.93 





24 





25 


BAL 


Emp 


„4.9 
"8.4 


0.96 


1 


37 





39 


non-BAL 


Obs 


"/it 


0.99 


9 


17 





22 


BAL 


Obs 
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Note. - 


See section 4 for 


a detailed explanation 


of the entries in 




this table. 

















a In this simulation the probability never reaches above 0.9, and so 
these values give the maximum range in 9 where the probability is over 
70% 



TABLE 3 

Simulation Results using R-S Tests 
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a 8A 
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Obs 


"8.4 
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0.99 


14 
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21 
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non-BAL 
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0.98 


13 
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21 
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a/it 


0.97 


15 
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27 


non-BAL 
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1.00 


7 


14 





21 
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a/it 


0.99 


17 


28 





29 


Note. - 
this table. 


See section 4 for 


a detailed explanation of the entries in 
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Fig. 1. — (Top) Observational data to constrain the a-8 relationship, with the derived linear fit. Here a is measured between 15 and 8.4 
GHz due to constraints on available data, and we assume that this can be extended to lower frequencies. Due to a lack of data and the fact 
that Doppler boosting does not likely effect the relationship above 30° , we use a fiat line there with some intrinsic scatter. (Bottom) Our 
monte-carlo simulation of 10,000 randomly distributed radio jets and their corresponding spectral index, determined by the observational 
data at the top. 




Fig. 2. — (Top) Simulation data from Wilman et al. (2008) showing the a-8 relationship for all FRI and FRII sources between 1.5 < z < 3.5 
with 51.4 greater than 10 mjy. Here alpha is measured between 151 and 610 MHz, due to issues with the Wilman et al. (2008) simulations 
at higher frequency. The linear fit is to the points between 10° and 20° before the data bottom out at a = —0.75 (a lower limit used in 
the models) and flatten out at low 8. The curved line is a linear fit to log(a + 0.75) and 8 between 5° and 20°. Below 10° the curved fit is 
used, above the linear fit is used. (Bottom) Our Monte-Carlo simulation of 10,000 randomly distributed radio jets and their corresponding 
spectral indicies, determined by the simulation data at the top. 



Fig. 3. — Comparison (via K-S tests) of simulations to observations of non-BALs (left) and BALs (right) for the observed ajf'4, using 
the observationally constrained model. The z-axis is the probability of the simulations matching the observations, as a function of 8 m i n 
(x-axis) and 6 m ax (y-axis). 




Fig. 4. — Comparison (via K-S tests) of simulations to observations of non-BALs (left) and BALs (right) for ct| |, using the semi-empirical 
model. The axes are the same as in Figure [3] 




Fig. 5. — Comparison (via K-S tests) of simulations to observations of non-BALs (left) and BALs (right) for auf, using the observational 
model. The axes are the same as in Figure [3] 
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Fig. 6. — Comparison (via K-S tests) of simulations to observations of non-BALs (left) and BALs (right) for a fit, using the semi-empirical 
model. The axes are the same as in Figure [3] 
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